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Abstract 

A simple criterion to optimise coarse-grainings for exact renormalisation group equa- 
tions is given. It is aimed at improving the convergence of approximate solutions 
of flow equations. The optimisation criterion is generic, as it refers only to the 
coarse-grained propagator at vanishing field. In physical terms, it is understood 
as an optimisation condition for amplitude expansions. Alternatively, it can be in- 
terpreted as the requirement to move poles of threshold functions away from the 
physical region. The link to expansions in field amplitudes is discussed as well. Op- 
timal parameters are given explicitly for a variety of different coarse-grainings. As a 
by-product it is found that the sharp cut-off regulator does not belong to the class 
of such optimal coarse-grainings, which explains the poor convergence of amplitude 
expansions based on it. 
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Introduction 

Flow equations |]J or exact renormalisation group equations are a promising tool for 
non-perturbative problems within quantum field theories. Based on the idea of integrating- 
out large momentum modes within a path integral formulation of quantum field theory, they 
permit to obtain an effective (or 'coarse-grained') theory for the remaining light degrees of 
freedom. Flows interpolate between some initial (classical) action, and the full (quantum) 
effective action. The particular strength of the formalism is its flexibility, allowing for 
systematic approximations without being tied to the small coupling region 

On the conceptual side, important progress has been made within the recent past, 
in particular through the consistent extension to (non-Abelian) gauge theories ||, which 
may play an important role in future applications to QCD. In turn, comparatively little 
attention has been paid to the interplay between approximate solutions to flow equations 
on the one side and the coarse-graining procedure on the other. Any explicit application 
of the formalism requires some approximation, simply because it seems very difficult to 
simultaneously solve infinitely many coupled partial differential equations. Approximate 
solutions have a finite domain of validity |4j], they may show a spurious dependence on the 
coarse-graining parameters at second order phase transitions || or first order ones ||[7|- 
and the corresponding expansions have a finite (often unknown) radius of convergence 
||. It has also been argued that the spurious scheme dependence can be used to assess 
the reliability of a given approximation ||, resulting in error bars describing the scheme 
and truncation dependence 0, or to obtain a better convergence employing a minimum 
sensitivity condition |J. 

The cross- dependence between the coarse-graining and a truncation of the flow has a 
very simple origin. Any coarse-graining couples to all operators in the effective action, 
with a strength depending on the details of the specific implementation. The sharp cutoff, 
for example, eliminates all momentum modes above some coarse-graining scale k, while 
smooth cutoff functions, like exponential, algebraic or mass-like ones, smoothly cut off 
the higher modes within some finite momentum interval about k. Hence, a change of 
the coarse-graining scheme alters the effective coupling to the operators in the effective 
action. A strong scheme dependence indicates whether some relevant operators have been 
neglected within a given approximation 0. Along similar lines it has also been argued 
that the smoothness of the coarse-graining can influence the contributions of higher order 
derivative operators [|]. 

More generally, this interplay suggests that a suitably chosen coarse-graining can sig- 
nificantly improve the convergence properties or the domain of validity of approximate 
solutions, which would be of great value for both analytical and numerical applications. In 
the present note we take profit of this freedom to propose a simple yet efficient optimisation 
criterion for Wilsonian flows, with the following characteristics: it optimises approximate 
solutions to flow equations; it is compatible with systematic approximation schemes like 
the derivative expansion; it only relies on a basic ingredient to the flow equation (the ef- 
fective coarse-grained propagator at vanishing field); it is generically applicable, i.e. not 
based on a specific theory. 
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Flows, coarse-grainings and propagators 

To begin with, we briefly review the basic ingredients needed for a Wilsonian flow equa- 
tion. The modern way to implement an exact renormalisation group procedure amounts 
to add a regulator term ~ / d d q(fi(—q) Rk(q 2 )4>(q) (for bosonic fields) to the action [|]. The 
operator Rk(q 2 ), sometimes referred to as the regulator scheme (RS), desribes the coarse- 
graining and introduces a fiducial scale parameter k, which is ultimately interpreted as 
a coarse-graining scale. This scale parameter induces a scale dependence, which, when 
written for the scale-dependent effective action T^, results in the flow equation 



Here, <p denotes bosonic fields and t = Ink the logarithmic scale parameter. The right 
hand side of ([!]) contains the regulator function and the second functional derivative 
of the effective action with respect to the fields. The trace denotes a summation over all 
indices and integration over all momenta. 

All specifications regarding the coarse-graining are given through the operator Rk(q 2 ). 
This operator can be chosen at will though within some basic restrictions which are briefly 
reviewed. First of all, it is required that Rk{q 2 — > 0) > 0. This ensures that the effective 
propagator at vanishing field Ak(q 2 ) = l/[q 2 + Rk{q 2 )} remains finite in the infrared limit 
q 2 — > 0, and no infrared divergences are encountered in the presence of massless modes. 
The second requirement is the vanishing of Rk in the infrared, Rk{q 2 ) — > for k — > 0. 
This guarantees that the coarse-grained generating functional of 1PI Green functions 
reduces to the usual generating functional of 1PI Green functions T = lim/.-^oIV The 
third condition to be met is that Rk{q 2 ) diverges in the UV limit k — > A. This way it is 
ensured that the microscopic action S = lim^A Tfc is approached in the ultraviolet limit 
k — > A. These conditions ensure that the flow (|l|) interpolates between the classical and 
the quantum effective action. For later convenience, we shall impose a further constraint 
on Rk, which, however, is nothing but a normalisation condition for the coarse-graining 
scale k, namely Rk{q 2 = k 2 ) = k 2 . If one would have chosen a regulator with Rt{q 2 = 
k 2 ) = ck 2 , then it is always possible to find a corresponding scale k c g = f(c)k such that 
R(q 2 = k 2 G ) = k 2 s . Hence, the last condition guarantees that (trivial) rescalings of k have 
been factored-out, and comparing different effective propagators is now sensible. 

It is important to realise that the integrand of the flow equation (|I|), as a function of 
momenta q, is peaked about q 2 ~ k 2 , and suppressed for large momenta. Consequently, 
at each infinitesimal integration step k — > k — Ak only a narrow window of momentum 
modes contribute to the change of — > r^Afc- Most importantly, modes with momenta 
g > no longer contribute to the running at the scale k. It is this property which justifies 
the interpretation of Tk as a coarse-grained effective action with modes q ^> k already 
integrated out. 

Let us now turn to the effective coarse-grained propagator Ak(q 2 ), or, for convenience, 
to its dimensionless inverse given by 




(1) 




(2) 
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and y = q 2 /k 2 . We find it convenient to write the regulator function Rk(q 2 ) in terms of a 
dimensionless function r(q 2 /k 2 ) as 

R k (q 2 )=q 2 r(q 2 /k 2 ) , (3) 

normalised as r(l) = 1. As a consequence of the above conditions, we can establish the 
following properties of P 2 (see Fig. 1). In contrast to the 'free' inverse propagator q 2 , the 
coarse-grained one is strictly positive for all momenta , P 2 > 0. This is due to the regulator 
term. For large momenta , P 2 grows linearly with y = q 2 /k 2 . 



Class 


P 2 (q 2 = 0) 


min 9 2 P 2 at 


la: 
lb: 
II: 


finite 
finite 
infinite 


q 2 = 



Tab. 1: Classification of different regulator schemes. 

For small momenta, r(y) diverges at least as 1/y, hence P 2 either approaches a constant 
> 0, or it diverges. We refer to the first class of regulators as mass-like, or Class I regu- 
lators, and to the second one as Class II regulators (see Tab. 1). Notice that their proper 
normalisation implies P 2 (q 2 = k 2 ) = 2 for any regulator. 
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Fig. 1: The effective inverse propagator for the different regulator classes la, lb and II 
(thick full lines). The physical normalisation implies that P 2 (y = 1) = 2 for all regulators. 
The dashed line corresponds to r = 0. 
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Optimisation 



The effective coarse-grained inverse propagator P 2 depends on the particular regular- 
isation scheme chosen. We can formulate a simple optimisation criterion, which is the 
requirement that the minimum of P 2 with respect to momenta be maximal with respect 
to the coarse-graining, 

C (mt = max ( minP 2 (y) ] = max ( min k 2 Ak(q 2 ) I . (4) 
p rs \ y >o yy, J rs \ q *>o\- vy J i J KJ 

Here, we denote those coarse-grainings as 'optimal' for which the maximum is attained. 
It is interesting to note that this criterion is based only on the effective propagator at 
vanishing field, that is, on the basic ingredient to the flow equation (p. Stated differently, 
no reference is made to a specific model or theory considered. This is desirable insofar as 
it renders the optimisation condition universally applicable. 

Some few observations concerning C opt can be made at this stage. For an optimal 
scheme, we have C opt = 2. This follows from combining r(l) = 1 with (|j), which implies 
mm y P 2 (y) < P 2 (y = 1) = 2 for any scheme. This also implies that the optimal P 2 reaches 
its (local) minimum at y — 1 and r = 1 for any optimised scheme. Such regulators are 
either of Class lb or Class II. A Class la regulator can only be optimal in the event where 
P 2 {y) has two degenerate minima, one at y = (by definition), and the other one at y — 1. 
Hence, an optimal Class la regulator is also of Class lb. The optimal value C opt = 2 cannot 
be reached for regulators, which are strictly Class la (P 2 {y = 0) is the global minimum). 

Threshold functions 

The optimisation criterion can be interpreted in more physical terms when generic 
threshold functions are considered. To give a particular example, consider a iV-component 
real scalar field theory in d dimensions. To leading order in the derivative expansion, 
we approximate the effective action as Tk = J d d x[\d^ L <j) a d tl 4> a + Uk(<f>) + 0(<9 4 )] and a = 
1, . . . , N. It is useful to introduce rescaled (dimensionless) variables as u (p) = U k /k d and 
p = \<p a ip a k 2 ~ d . The flow equation for the first derivative of the effective potential v! is a 
second-order partial differential equation, to wit 

d t u' = -2v! + (d- 2)pu" - 2v d (N - l)u" ^(wi) - 2v d {Zu" + 2pu"')£f(uj 2 ) (5) 

with v J 1 = 2 a!+1 7r a! / 2 r[|]. Here, the relevant dimensionless amplitudes are uj\ = vl and 
u>2 = u' + 2pu", and the threshold function £f(u>) reads 

dyr'(y)y 1+d ' 2 — j- (6) 

yP + UJ) Z 



Flow equations of the form (|5|) have been intensively studied in the literature |T0|-|T2, 

In the general case, an Ansatz for r& need not be restricted to the leading order in the 
derivative expansion. At higher order, one finds coupled sets of partial differential equations 
for the coefficient functions of the higher derivative operators, involving more complicated 
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threshold functions. They are obtained from the flow equation (|lj), or functional derivatives 
thereof, within specific models and approximations. A generic threshold function has the 
following structure, 

dyK ni _ ni (y) [P\y) + . . . [P\y) + (7) 

and all rii > 0. Typically, the amplitudes Wi describe (field-dependent) mass terms uj ~ 
ml(<f))/k 2 , which grow large once the coarse-graining scale falls below the mass threshold, 
k <C mfc(0). This leads to the decoupling of the corresponding modes (hence the name 
'threshold' functions). In dimensionful units the physical amplitudes are Q = k 2 u ~ m 2 (<j)). 
The explicit form of the kernel K ni) ___ tni depends on the specific quantity under investigation. 
In the example above, the threshold function (|6|) is of the form (|7]) for the specific kernel 
K 2 (y) = -r'(y)y 1 + d / 2 . 

As an immediate consequence of the properties of the basic flow equation ([!]) we can 
conclude that the kernel of (0) is well-defined, finite and peaked as a function of momenta. 
Given the properties of P 2 , we deduce that such threshold functions always have a pole 
on the negative u;-axis, located at uj po \ c = —mm y > P 2 (y) < 0. The location of the pole 
depends on the particular coarse-graining chosen. Notice that the pole cannot be shifted 
arbitrarily far away towards the negative axis. Requiring that the pole of any threshold 
function be as negative as possible corresponds precisely to the condition (^). Hence, the 
maximally attainable point is given by C opt = —u) po \ e = 2. The condition (|j) applies to 
any threshold function of the form (0), and therefore holds for any amplitude uj > uo po \ e . In 
dimensionful units, the pole is located at fl po i e = k 2 u) po \ e < 0. For k close to some ultraviolet 
scale A the pole is far away on the negative axis and outside the physical domain. It is 
moving inward for decreasing k, ultimately reaching Q po \ c = in the infrared limit, where 
the potential becomes convex. 

The pole structure of threshold functions is closely linked to the flattening of the non- 



convex parts of effective potentials in the limit k — > 0, like Uk in the example above |13 
The convexity is a consequence of the 1PI generating functional being obtained by a Legen- 
dre transformation. For finite k > 0, though, convexity is not required. In the non-convex 
region, the amplitudes uj\ ~ A;~ 2 _1 <9£4/<90 or u> 2 ~ k~ 2 d 2 Uk/d(j) 2 turn negative. In par- 
ticular, for k — > 0, they approach the poles of threshold functions. Hence, choosing a 
coarse-graining for which this pole is the most negative corresponds to choosing a coarse- 
graining for which the transition towards a convex effective action is the smoothest. Such 
a choice is particularly useful for the stability and an improved convergence of numerical 
implementations of flow equations. 



Amplitude expansions 

An alternative physical explanation of the optimisation condition is given in terms of 
an expansion in (small) amplitudes. With amplitude expansion we mean any expansion 
in the amplitudes to in units of P 2 . Such expansions always exist, simply because P 2 
is strictly non-zero. Then, all information about the regulator scheme is encoded in the 
related expansion coefficients, introduced in 0. Here, they can be written as 
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a n = / dyK[r]P- n (y) 
Jo 



The explicit form of the kernel K depends on the specific quantity studied. For the ex- 
ample given above, the threshold function is Taylor-expanded in powers of the amplitudes 
£f(uj) = J2m=i if a 2m-d (—uj)™ 1 ^ 1 where the coefficients a n are those defined in (||) with the 
particular kernel 

KM = -- ^ (9) 

L J 2[l + r(y)] 1 +<V 2 ' { ) 

normalised such that a = 1 ||. 

The main characteristics of amplitude expansions are deduced from the knowledge 
that the kernel K for general expansion coefficients (|8|) is finite, peaked, and suppressed 
for sufficiently large momenta q 2 /k 2 . Conceptually speaking, it would be most desirable if 
only a few terms in such a series have to be taken into account, and higher order corrections 
be small. Hence, we look for coarse-grainings for which the main physical information is 
already encoded within a few leading order terms. For fixed amplitudes, this amounts to 
the requirement that the higher order coefficients a n be as small as possible, or, that the 
radius of convergence for amplitude expansions be maximal. The (dimensionless) radius 
of convergence C is defined as the ratio a n /a n+ 2 of two successive expansion coefficients in 
the limit n — > oo, 

C = lim . (10) 

This limit is computed as follows. First consider the definition (||D and notice that the 
term P~ n will strongly suppress a n in the limit n — > oo because P 2 is strictly positive 
and diverging for large momenta (see Fig. 1). The sole contribution to the integrand will 
then come from the minimum of P 2 where the integrand is the least suppressed. Second, 
taking into account that the kernel K is well-behaved, we can conclude that the radius of 
convergence is given by the 'saddle point approximation' to fllPf) , namely 

C = mm P 2 (y) . (11) 

The radius of convergence does not depend on the kernel K, but still on the parameters of 
the coarse-graining. The optimisation criterion for amplitude expansions C op t = max^s C 
coincides with the optimisation condition introduced earlier. The finiteness of C is due 
to the pole structure of threshold functions. This implies that the radius of convergence 
cannot be increased by changing to another expansion point. The optimal expansion point 
is to — 0. 

The optimal coarse-grainings have another interesting interpretation. Consider again 
the definition @ for the expansion coefficients, and let us split them as a n = a< + a> into 
its low momentum contributions a< = Jq 1 dyK [r]P~ n and its high momentum contributions 
a n = Ji°° dyK[r]P~ n . We assume that the coarse-graining scheme is parametrised by some 
parameter b. Let us define implicitly the parameters b n from the requirement that the hard 
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modes and soft modes give equal contributions to the coefficients a n , hence < This 

is certainly possible for n sufficiently large. The optimal choice for the coarse-graining 6 opt , 
as defined by the optimisation condition, can be shown to correspond precisely to the limit 



Field amplitude expansions 

It is worth pointing out a crucial difference between expansions based on the amplitudes 
u>i, as considered here, and expansions based directly on field amplitudes <p a or p (sometimes 
also denoted as a local polynomial approximation), which have often been used in the 
literature 



We consider the iV-component scalar field theory with the flow equation (|^) as example. 
Expansions in field amplitudes p of the effective potential read u' = Y.n=i w.^ n (P ~ ^o)™? 
where Ao 7^ corresponds to the chosen expansion point The partial differential equation 
(|) reduces to a set of coupled ordinary differential equations (3 n = d\ n /dt for the (M + 1) 
couplings A n . These /3-functions depend on the coarse-graining only through the threshold 
function if (u) and derivatives thereof, evaluated at the expansion point uj(X ). The relevant 
amplitudes are uj\ = u' and oj 2 = u' + 2pu", and the optimal expansion point is to\ — u 2 — 0. 
When expressed in terms of field amplitudes p, these expansion points are not equivalent .[] 
An expansion in p therefore has a smaller radius of convergence than expansions in the 
amplitudes u>i. 

Interestingly, there are two exceptions in the present example. The first one is the limit 
iV — > 00, in which case only the amplitude u>i = u' remains. The expansion about the local 
minimum of the scaling potential p\ with uJi(pi) = is the optimal choice and equivalent to 
the expansion in the amplitude iO\. A second exceptional case concerns N = 1, where only 
the amplitude u 2 = u' + 2pu" appears, suggesting that an expansion about the inflection 
point p 2 with u 2 (p2) = has best convergence properties, as long as u> is Taylor-expandable 
about p2- 

It was already argued that expansions about the local potential minimum p\ show bet- 
ter convergence than expansions about vanishing field Ao = |14|| . The present discussion 
clarifies from a more general perspective why an appropriately chosen field expansion point 
does improve the convergence of the entire series. While the optimal expansion point for 
amplitude expansions is always Ui = 0, the optimal expansion point for field amplitudes 
Popt depends on the particular theory studied (in the example discussed here on the number 
of real scalar fields N). Clearly, all these expansions, either in p or in the amplitudes u>i, 
are optimised through an appropriate choice of the coarse-graining, as defined in (|). 



1"For expansions about p = we write v! = J2n=o ^\^nP n instead. 

^The condition u)\ = u) 2 = over-constrains the choice for an expansion point Ao as it implies 
either v! = u" = 0, or u'(0) = and u"(0) = finite. 



7 



Applications 

We now turn our attention to a number of specific coarse-graining functions. Our 
list of regulators is by no means exhaustive, though it covers the main regulators used 
and discussed in the literature. It includes exponential, algebraic and mass-like ones, 
and the sharp cut-off for comparison (see Tab. 2; the parameter c ensures the consistent 
normalisation of the scale parameter k). They all differ in the way how high momentum 
modes are cut off, with the parameter b measuring the relative 'smoothness' - the smaller 
b the 'smoother' the coarse-graining. 



Regulator 



Parameter range 



Class 


b pt 


Copt 


II* 


1.44 


2 


I 


3.92 


2 


I 


1.92 


2 


II 


2 


2 


I 


4.40 


1.65 


II* 


2 


2 


II 




1 



'"exp 


exp cy b — 


1 -1 

1 




'"mexp 


b [exp cy - 






'"mod 


exp 


l(y + (b-l)y b )- 


'"mix — exp 


-Uv a -v' a \ 




'"step 


26-2 
b 


yb-2 


expcy b 1 




'"power 


y-b 








'"sharp 


1 








e(y- 


1) 1 







c = In 2 
c = ln(l + 

c = In 2 

a > 

c = ln^ 





1 < b < oo 
1 < b < oo 



Tab. 2: The definition of various classes of regulator functions and their optimal parameter 
choices, including the sharp cutoff for comparison. Class II* denotes Class II regulators 
except for b = 1, where they correspond to Class I. 



More specifically, the regulators r exp , r mod , r step , r mexp for b ^ and r mix for a ^ 1 are 
exponential ones, that is, exponentially suppressed for large momenta. In contrast, r power , 
fmexp for 6 = 0, and r mix for a = are algebraic regulators. In the classifaction of Tab. 1, 
the mass-like exponential regulator r mexp , the modified exponential regulator r mo( j and r step 
are Class I (mass-like) regulators, as are the exponential one r exp and the algebraic one 
r power for b = 1. On the other hand, the mixed exponential regulator r mix , the sharp cut-off 
'"sharp, and r exp and r power for b > 1, are Class II regulators. 

Some cross-dependences amongst these functions are worth noticing. For b — > oo, 
'"exp, '"mexp, '"mod and r power approach the sharp cut-off limit r s h arp , while r step approaches 
~ y~ l 6(l — y), which corresponds to a mass term which is cut-off for momenta above k, 
-Rstep, k ~ k 2 9(k 2 — q 2 ). The two-parameter family of mixed exponential regulators r mix , 
which obey r mix (l/?/) = l/r mix (?/), contain the algebraic regulator r power as the limiting 
case a — > 0. r exp and r mo d for b = 1 coincide with r step for b = 2. Also, r mexp for b = 0, and 
r step for b = 1 correspond to r power for b = 1. 

We now distinguish the different coarse-grainings by their associated radii of conver- 
gence. In Fig. 2, we have displayed C, as defined in (0), for different coarse-graining 
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schemes and as functions of the coarse-graining parameter b. In the sharp cutoff limit 
b — > oo the value C s har P = 1 is attained. For the classes of regulators considered here, we 
notice that the maximal radius is reached for specific smooth coarse-grainings. The radii 
vary between 1 < C < 2 (though it is possible to find regulators such that C is arbitrarily 
small). Optimal parameters are found for both Class I and Class II regulators (see Fig. 3). 
The most efficient way for determining optimal parameters for all Class lb and Class II 
regulators is the following. We have established that P 2 {y) takes its maximum at y = 1 
and r = 1 for an optimal coarse-graining. The vanishing of dP 2 / dy = 1 + r + yr' at the 
extremum implies r' t (l) = —2. This allows to fix the optimal coarse-graining parameters. 




Fig. 2: The dimensionless radii of convergence for amplitude expansions (a.k.a. location 
of the pole for threshold functions) as a function of the coarse-graining parameter for the 
regulators r power , r exp and r mo d (left panel) and r m i x with a = 0, |, 1 and 2 (right panel). 



For r exp one finds 6 exp = l/ln2, for r power the radius of convergence reads C — 6(1 — 
b) l / b ^ 1 which attains its maximum at 6 pow er = 2. The optimal parameter for r m j x is 6 m i x = 2 
for all a. The value for 6 mexp is obtained numerically by solving 2b = (6+1) ln(6 + 1), and 



the one for 6 mod reads 6 mo d = [1 + 5^5] [| + \J\ — 1/(1 + j^) 2 ] (see Tab. 2). The regulator 
r step does not attain the maximal radius C = 2. The reason behind it is that Ps tep (y) has 
two local minima for b ~ b opt , one at y ~ 1, and another one at y = 0. The minima 
are degenerate for b = b opt , and their competition only allows for a sub-optimal value 
Cstep — 1-65. 

The fact that C s har P = 1 is half as big as the possible optimal value explains why 
amplitude expansions of sharp cut-off flows have a notoriously bad convergence behaviour 
[TT] | . This property is not a particularity of the sharp cut-off; the smooth cut-offs r powcr , 
r exp and r mexp (all for b = 1) also have a radius of convergence C — 1. A significant 
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improvement is achieved only for the optimal coarse-grainings, which all belong to smooth 
cut-offs. 




1 - 

o L^l — i 1 1 1 1 

0.5 1 1.5 2 2.5 3 

q 2 /k 2 

Fig. 3: The optimised effective inverse propagators for different regulators (thick full lines). 
From the left-lowest curve clockwise: r mexp , r mod , r exp , r power , r mix (a = \), r mix (a = 1) 
and r m j x (a = 2). The dashed line corresponds to r = 0. 



The present discussion can be extended in several directions. Fermionic fields can be in- 
cluded using a regulator term like ~ J d d q%[){—q) Rf^q) ip(<l), with R F k (q) = q 2 rp(q 2 /k 2 ). 
Effectively, it appears in threshold functions in the combination Pp = y[l + r-p(y)} 2 , which 
is the fermionic analogue of (|2|). Hence, the optimisation condition for fermionic degrees 
of freedom is the same as (^), replacing P 2 by P|, yielding CF, op t = 4. For gauge fields, 
the flow equation is amended by a (modified) Ward or a (modified) BRST Identity, which 
ensures that physical Green functions obey the usual Ward or BRST Identities 0j. The 
criterion (§) is compatible with such an additional constraint. Also, wave function renor- 
malisations can be taken into account in the usual manner, setting Rf, = Zkq 2 r(q 2 /k 2 ). 
This modification will not change the optimisation as discussed in the present note, and 
allows for a systematic inclusion of higher derivative operators. The same holds for field 



theories at finite temperature within the imaginary or the real-time formalism flT5] . 

In conclusion, we have presented a simple and generic optimisation criterion for exact 
renormalisation group equations. It is interesting to mention its close link to a minimisa- 
tion of the spurious scheme dependence of flow equations, a more detailed discussion of 
which shall be presented elsewhere [fTEfl . 

The author thanks Ulrich Ellwanger, Filipe Freire, Jan M. Pawlowski, Janos Polonyi 
and Christof Wetterich for useful discussions. 
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